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Studying the energy partioning published by Q we show that the presented results do not fulfill 
the sum rule due to energy conservation. The observed fluctuations of the energy conservation test 
point to a numerical problem. Moreover, analysis of the binding energies show that the fragment 
recognition algorithm adopted by Campi et al. leads with a sizeable probability to fragments 
containing up to the total mass even for excitation energies as large as 3/4 of the total binding. 
This surprising result points to another problem since the published inter-fragment energy is not 
zero while a unique fragment is present. This problem may be due to either the fragment recognition 
algorithm or to the definition of the inter and intra-fragment energy. These numerical inconsistencies 
should be settled before any conclusion on the physics can be drawn. 

In a recent article^ X. Campi et al. present results from numerical experiments on the liquid gas phase transition 
of a Van der Vaals fluid using a molecular dynamics simulations. They acknowledge the fact that in a microcanonical 
ensemble the kinetic energy fluctuation (er^-/T 2 ) can be used to reconstruct the heat capacities in the context of phase 
transitions as proposed in ref. [2], but they criticize the experimental evaluation of such a quantity from the observed 
partitions as used in ref. 0. 

Before discussing the actual criticisms of Campi et al. to our work, we would like to point out that the numerical 
results presented in ref. are inconsistent. 
Campi et al. present 4 quantities, 

• the exact kinetic energy K = E — V, where E is the total energy and V the total interaction energy, 

• the quantity supposed to be the approximation to the kinetic energy introduced in ref. K 1 = E — Bi 
where the sum runs over the fragments and where Bi are the fragment ground state energy , 

• the inter-fragment interaction Vinter = 53j<j Vij 

• and the intra-fragment energy Vintra = Ylj Vj which they show as a configurational excitation energy AV = 

vr, n, 

Since V = Vinter + Vintra the above quantities should at any time fulfill the exact sum rule 

S = K - K' + V mter + AV = (1) 

The right side of figure 1 shows the sum rule as obtained from the data of ref. l] , which are shown in the left 
part of the same figure. 

Fig 1 shows that in the Campi et al. calculations the sum rule S is not zero, suggesting a possible numerical 
mistake. If the error on the average < S >— —8.871 is of the order of a few percent, the more worrying, considering 
that the aim of the Campi et al. paper is to discuss fluctuations, are the enormous fluctuations presented by this 
quantity which should be exactly conserved in time. Indeed the observed variance of S is 3819 while the one of K' 
criticized by Campi et al. to be too large is only 784. The problem can be more easily understood looking at figure 2 
which is the same as figure 3 of ref.pj on an enlarged time scale. In this figure we can see that the fluctuations of AV 
and Vi n ter are approximately in phase. If we consider that the fluctuations of K are very small as shown in figure 1, 
for the sum rule to be exactly fulfilled, K' should be positively correlated with AV and Vinter which is clearly not 
the case in the calculation of Campi et al. 
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FIG. 1: Left side: Fig 3 of [l| showing from bottom to top Vinter , AV , K' shifted by 100 units , K shifted by 250 units. Right 
side: the corresponding energy sum rule S — K — K' + Vi n t er + AV. 




FIG. 2: Same as the left part of Figure 1 on an enlarged scale, showing from bottom to top Vinter, K' shifted by —100 units 
and AV. 

Therefore, before entering into any discussions about the physical meaning of Campi et al. findings, one should be 
sure that the numerics is under control. 

In order to understand better what might be going on, we have used K, Vi n t er and AV combined with the Lennard 
Jones mass table^j) and the excitation energy reported in pj to disentangle the various components entering in the 
calculation and in particular the total binding energy of the partitioned system, B = ^ i Bi. We have thus used for 
the total energy E = E* + B$4 = —73.62 since E*/N — 4 according to and since £?64 = —329.62 according to 
the mass table|4(. Then we can get V from K = E — V and thus V intra from V = V inter + V intra and finally B from 
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FIG. 3: Left side: from bottom to top, the ground state binding energy of the N=64 system B@4 from ref.0|, the binding energy 
of the partitions B, the intra- fragment energy Vintra and the inter- fragment energy Vinter, for the same data as in figure 1. 
Right side: same as the left side, but on an enlarged time scale. 
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Ay = Vintra — B. The results are shown in figure 3. Before discussing the fact that B is a very poor approximation 
of Vi n tra m the calculations performed by Campi et al., one is surprised by the fact that B is reaching exactly 
£?64 = —329.62 at time step 711 in the figure 3 of ref. pj. This means that there is a non negligible probability that 
the system is composed by a unique fragment containing exactly all the particles (for comparison = —323.49). 
For this peculiar single-fragment event Campi et al. report a non zero inter-fragment energy (Vinter — ~ 29), which 
again points to a numerical problem. This might be related to the sum rule problem discussed above but it is likely 
to be a different one, since the first problem is related to K' while in the second one the information about K' is not 
used. One of the problems may lie in the fragment definition algorithm employed by Campi et al. In this respect it 
may be worthwhile to notice that, even if the classical Lennard Jones molecular dynamics is in itself an exact well 
controlled model, on the other hand the definition of physical clusters in these kind of models and their relationship 
with measured fragments is subject to many controversies 

In conclusion, our study of the results published by Campi et al. clearly show two numerical problems: 

• the sum rule taking care of the energy conservation in the energy partition is violated and presents fluctuations, 

• the analysis of the total binding energy of the partitioned system shows that even at excitation energies above 
the critical point, the fragment recognition algorithm used in ref j| predicts the existence of residues of very 
large masses, containing up to the 64 particles of the considered system, and pointing to a problem in the 
definition of the inter-fragment energy since in the case of a unique fragment this latter is not zero. 

These numerical inconsistencies should be settled before any conclusion on the physics can be drawn. 

To progress on this issue, we have repeated the Campi et alanalysisjlj of the independent fragment approximation 
for the evaluation of the interaction energy on another well controlled exact model. We show in a forthcoming paper 
that this approximation is accurate provided that the required consistency tests |6( are taken into account. 
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